Pseudogap of Color Superconductivity in Heated Quark Matter 
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We show that the pseudogap of the quark density of states is formed in hot quark matter as a 
precursory phenomenon of the color superconductivity on the basis of a low-energy effective theory. 
We clarify that the decaying process of quarks near Fermi surface to a hole and the diquark soft mode 
(qq)soft is responsible for the formation of the pseudogap. Our result suggests that the pseudogap 
is a universal phenomenon in strong coupling superconductors. 
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It is an intriguing problem to examine how rich is the 
phase structure of the high-density (p) QCD (Quantum 
Chromodynamics) matter at vanishing or moderate tem- 
perature (T). It is now believed that quark matter at an 
extremely high density undergoes a Cooper instability 
leading to the color superconductivity (CS) 0, S IE l3 ■ 
As might be ensured by the asymptotic-free nature of 
QCD, the quark matter at extremely high densities is 
usually treated as a Fermi liquid, and the mean-field 
(MF) theoretical approach a la BCS theory is einployed; 
see, for reviews on recent exciting development, |3- 

Various complications are, however, expected on the 
nature of quark matter at intermediate baryon densities 
or chemical potential /z where the strong coupling nature 
of QCD may show up and invalidate the MF approxi- 
mation (2| ■ The strong coupling may make the so-called 
Ginzburg region so wide that precursory fluctuations of 
the quark-pair field can have a prominent strength and 
may give rise to physicallv significant effects above the 
critical temperature (Tc)jg]. In this Letter, we shall show 
that the fluctuation effects of the CS can be so signifi- 
cant that such quark matter may share some basic prop- 
erties with the cuprates of the high- Tc superconductivity 
(HTSC), having a non- Fermi liquid nature above Tc. 

One of the typical non-Fermi liquid properties of the 
cuprates is the existence of the pseudogap, i.e., an anoma- 
lous depression of density of states (DOS) N{lo) as a func- 
tion of the Fermion energy lo around the Fermi surface 
above Tc. Although the mechanism of the pseudogap as- 
sociated with the HTSC is still controversial, precursory 
fluctuations of the pair fleld and the quasi-two dimen- 
sionality of the system seem to be basic ingredients to 
realize the pseudogap Q- In this work, we shall show 
that the pseudogap of the quark density of states exists 
as a precursory phenomenon of the CS in a considerable 
range of T above Tc even in the three-dimensional sys- 
tem. We calculate the quark self-energy incorporating 
the pre-forming pair-field of the CS for the first time in 
the T-matrix approximation (T-approximationU0,^3|)j 
where the amplitude fluctuations are assumed to be dom- 
inant over the phase fluctuations of the pair fleld 0, [ll| ■ 
Our result shows that the pseudogap may appear also in 
a relativistic system solely owing to the amplitude fluc- 



tuations, and suggests that the pseudogap phenomenon 
is universal in strong coupling superconductors, irrespec- 
tive of the dimensionality as suggested from the study 
of nuclear matter 0. Our results should also provide an 
insight into the physics of proto-neutron stars and heavy- 
ion collisions as well, where the p is relatively low and the 
effect of finite T plays an important role. 

To describe a system at relatively low T and p, it 
is appro pria te to adopt a low-energy effective theory of 
QCD m Il2l Il3| . Here we employ a simplified version of 
the instanton-induced interaction in the two-flavor case 
known as the Nambu- Jona-Lasinio modeliflHi with the 
scalar-diquark interaction in the chiral limit, 



(1) 



where = Cip'^{x), with C = 17270 being the 

charge conjugation operator. Here, T2 and mean the 
antisymmetric flavor SU(2) and color SU(3) matrices, re- 
spectively. The coupling Gs and the three dimensional 
momentum cutoff A = 650 MeV are determined so as 
to reproduce the pion decay constant = 93 MeV and 
the chiral condensate {ip'ip) — (— 250MeV)^ in the chi- 
ral hmit 0. We choose Gc so as to reproduce the 
phase structure calculated using the instanton-induced 
interactionjll, i.e., Gc = S.llGeV-^ [3. We neglect 
the gluon degrees of freedom, especially their fluctuation, 
which is known to make the CSphase transition flrst or- 
der in the weak coupHng regionfalT^: notice that the CS 
is a type I in this regime. However, as is emphasized in 
Ref. [13, nothing deflnite is known on the characteris- 
tics of the CS in the intermediate density region. In this 
work, simply assuming that the fluctuations of the pair 
fleld dominates that of the gluon fleld as is the case for 
type II color superconductors, we examine the effects of 
the precursory fluctuations of the diquark pair fleld on 
the quark sector in the T-approximation. 

The DOS N{u!) is calculated from the spectral function 
A{k,uj) of a single quark, which is deflned through the 
spectral representation of the retarded Green function of 
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the quark field; 



r^' ■"■''•r'\ , (2) 

J UJ-UJ'+IT]' ^ ^ 



and accordingly, A{k,uj) — — I/tt • IniG-'^(A;, w) = — I/tt • 
(G^ — 7°G^^7'')/2i. From the rotational and parity in- 
variances, the spectral function has the following matrix 
structure: A{k,uj) ~ pQ{k,uo)^^ ~ p^{k,uj)k-^ + ps{k,u}), 
where k = k/\k\ and pa (a = 0,v,s) still have color 
and flavor indices. Since the quark number is given by 
N ^ J d^-x.{^^'^'ip), the DOS is solely given by pQ{k,uj), 



N{lu) = 4 



d^k 
(2^ 



Trc,f [po{k,uj)] , 



(3) 



with Trc,f denoting the trace over color and flavor indices. 

The G^ in Eq. Q is given by the analytic continua- 
tion of the imaginary-time (Matsubara) Green function 
Q, which obeys the following Dyson-Schwinger equation 

g{k,uj„) ^ go{k,u;n){l + ^{k,uJn)g{k,uJn)}, (4) 

where Go^k, ujn) and S(fc, uJn) denote the free Green func- 
tion and the self-energy in the imaginary time, respec- 
tively. In the normal phase, Go is reduced to Ga{k, LUn) — 
[{iujn + p)^^ — k ■ 7]~^ with the Matsubara frequency 
LJn = (2r7, -|- 1)'kT for fermions. 
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As was shown in [Gj, the fluctuating diquark pair field 
develops a collective mode (the soft mode of the CS) at 
T above but in the vicinity of Tc, in accordance with the 
Thouless criterion T^. Our point in this work is that the 
soft mode in turn contributes to the self-energy of the 
quark field, thereby it can modify the DOS so much to 
give rise to a pseudogap. The quark self-energy E owing 
to the soft mode may be obtained by the infinite series 
of the ring diagrams shown in Fig. 1; 



S(fc, uj„ 



d^ki ~ 



-mc{i + GcQ{k,v^)y 



(5) 
(6) 



with the lowest particle-particle correlation function 
Q{k,Vn) and Vn = 2rmT being the Matsubara fre- 
quency for bosons. Notice that the thin quark lines in 
Fig.l are the free Green function, so we have taken the 
so-called non-self-consistent approximation, on which we 
shall make a comment later. 

Inserting Eqs. lO and © into Eq. Q and performing 
the analytic continuation to the upper half of the com- 
plex energy plane, we obtain the retarded Green function, 
G^(fe,w) = (G(7^(fc,w + ir;) -S«(fc,cj))-i, with 



E«(p,/) 



dw -Im :i^{p + q, uj) Eq^^ - Q • 7 



dh 
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FIG. 1: The Feynman diagrams representing the quark Green 
function. The thin lines represent the free propagator Qo, 
while the bold ones the full propagator Q. 



where Eg = |q|, Y.^{k,Lu) = T,{k,uj„)\iu„=u+iri and 

, i^n)\iu„=uj+iri- The matrix structure 
of the self-energy is the same as the spectral func- 
tion, E^{k,uj) = So(A;,cj)70 - Y.^{k,uj)k ■ 7 + Es(fe,w) 
EE 70(I]_A_ + S+A+), where Azp = (1 ± 7% • k)/2 
denotes the projection operators onto the positive and 
negative energy states, and accordingly Ezp = Eq T ^^v 
represents the self-energies of the particles and anti- 
particles, respectively; notice that Es(fc,w) = in the 
chirally restored phase in the chiral limit. Then it 
can be shown that po is expressed as a sum of the 
positive- and negative-energy parts; pQ{k,uj) — — l/27r • 



E„^±ImE„/{i?„(fc,a;)2 + (ImE^)^}, where R±{k,uj) - 
uj±\k\+ p~ ReJ:±{k,uj). 

A remark is in order here: The Thouless criterion 
mentioned above tells us that the denominator of E^, 
1 -I- GcQ{O,0), vanishes dX T — Tc because of the 
self-consistency condition for the diquark condensate at 
T = TcQ, where Q{k,uj) = Q{k,Vn)\iv^=uj+irj- 

For the numerical calculation, we employ the following 
cutoff scheme [T3.IT9l |: First we notice that the imaginary 
part of Q is free from ultraviolet divergences. Therefore 
we first evaluate the imaginary part exactly without in- 
troducing a cutoff, and then use the dispersion relation 
with the imaginary part just obtained to calculate the 
real part introducing a three dimensional cutoff A. The 
imaginary part is nicely found to have the following com- 
pact form 



Im Q(fe,w) = -Nf{Nc ■ 
X 9{\uj + 2p\- k)lii 



{uj + 2py 



- 1)T 

' 27rfc 

cosh(w -f fc)/4T 
cosh(w - fc)/4T 
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2^1 + k) In 



1 



-(uj+k)/2T 



(8) 



We remark that each term has the respective interpre- 
tation in terms of the kinetic processes 0. We empha- 
size that the compact expression of Im Q{k, lu) above ex- 
tremely improves the efficiency of the numerical calcula- 
tions. Nevertheless the calculations still involve multiple 
integrations for obtaining and hence pQ. 

Since pQ{k,Lu) for to > —fi is well approximated solely 
by the positive-energy part, we first see the charac- 
teristic properties of E_. Figure |5] shows a typical 
behavior of the real and imaginary parts of S_ with 
k = kp at fj. = 400MeV and the reduced temperature 
e = (T - Tc)/Tc = 0.01 ; we remark that Tc = 40.04MeV 
in the present case|{|j. From the figure, one can see that 
ReS_ shows a rapid increase around the Fermi energy 
uj = 0. The quark dispersion relation u = LO-{k) there- 
fore also shows a similar behavior around the Fermi sur- 
face; uj-{k) is the solution of R^{k,uj) = 0. Hence the 
density of states proportional to {duj- /dk)~^ becomes 
smaller near the Fermi surface, which suggests the exis- 
tence of a pseudogap, provided that the imaginary part 
ImE_ is neglected, which will be discussed shortly. One 
can also see that a;_ [k = kp) ~ k—p, since Rel]_ at a; = 
is vanishingly small, which will be found important for 
the pseudogap formation around the Fermi surface. 

A numerical calculation shows that as the momentum 
k is varied from kp, the peak of |Im$]-| at w « seen in 
Fig|21 moves along to k ~k + jjL^^. This means that the 
quasiparticles with this energy are dumped modes. Fig- 
ure ^ tells us that ImE_ describes a decay process of a 
quark to a hole and a diquark, q^h-l-(qq), where the hole 
is on-shell with a free dispersion relation ujh = p, — \kh\. 
The essential point for the pseudogap formation is that 
the above process is enhanced when the diquark (qq) 
makes a collective mode, which we have emphasized is 
the case; the diquark soft mode (qq)soft has a promi- 
nent strength at vanishingly small energy and mo- 
mentum ks near T,.. Owing to the energy-momentum 
conservation, the energy-momentum of the decaying par- 
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fes ~ 0. It means that when the decaying 
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FIG. 2: The self-energy E- with k = kr sX p. = 400MeV 
and e = {T — Tc)/Tc = 0.01. One observes a peak in ImS_ 
at Lj = and a rapid increase of ReE_ at the same uo. 
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FIG. 3: The spectral function po ai ^ — 400MeV and 
e = 0.01. The upper figure is an enlargement of the one near 
the Fermi surface. The peaks a.t uj = k ~ p and lo = — — p 
correspond to the quark and anti-quark quasiparticles, respec- 
tively. Notice that there is a depression around = 0, which 
is responsible for the pseudogap formation. 



particle has almost the same energy as a free quark has, 
|ImI]_(fe,w)| has the largest value. 

The spectral function pQ{k,uj) is shown in Fig|31 at 
the same p. and e as those in Fig|21 One can see two 
families of peaks around lo = a;_(fc) k, k ~ p and 
Lj = —k — p, which correspond to the quasiparticle peaks 
of the quarks and anti-quarks, respectively. A notable 
point is that the quasiparticle peak has a clear depres- 
sion around lu — 0, i.e., the Fermi energy. The mech- 
anism for the depression is easily understood in terms 
of the characteristic properties of the self-energy men- 
tioned above: In fact, R-{kp,LU ~ /i — kp) ~ and 
|ImI]_(fe, 0^)1 becomes large when lu « —k + p. Thus 
Po{k ~ kp,uj ~ 0) ~ — l/(27rlml]_(fci?, w ~ 0)), which is 
found to be suppressed. 

Integrating po, one obtains the DOS N{lu): Figure^ 
shows the DOS at ^ = 400MeV and various values of 
the reduced temperature e together with that of the free 
quark system, Nq(lu). As anticipated, one can see a re- 
markable depression of N(lu), i.e., the pseudogap, around 
the Fermi energy lu = 0; N (lu) / No{ui)\^=o ~ 0.64 at 
e = 0.01. One sees that the smaller e, the more re- 
markable the rate of depression. The clear pseudogap 
structure survives even at e = 0.05. One may thus con- 
clude that there is a pseudogap region within the normal 
phase above up to T = (1.05 ~ l.l)Tc at p = 400MeV, 
for instance. A numerical calculation shows that e de- 
pendence of the width of the pseudogap region hardly 
changes for 320MeV< p < 500MeV. 

We notice that the pseudogap region obtained in the 
present work is more than one order of magnitude wider 
in the unit of e than in the nuclear matter '3| where the 
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FIG. 4: Density of states at jj. — 400MeV and various e = 
(T — Tc)/Tc- The Dotted line shows that of the free quarks. 
A clear pseudogap structure is seen, which survives up to 
£ 0.05. 



clear pseudogap is seen up to e w 0.0025. This is just 
a reflection of the strong coupling nature of the QCD 
at intermediate density region. Our result obtained for 
a three-dimensional system tells us that a considerable 
pseudogap can be formed without the help of the low- 
dimensionality as in the HTSC and that the pseudogap 
phenomena in general may be universal in any strong 
coupling superconductivity. 



In this Letter, we have found that the pseudogap can 
be formed as a precursory phenomenon of the CS in 
a rather wide region of T above Tc. This may imply 
that quark matter shares some basic properties with the 
cuprates of the HTSC. It should be noted that our work 
is the first calculation to show the formation of the pseu- 
dogap in the relativistic framework. 

In the present work, we have employed the non-self- 
consistent T-approximation, while the self-consistent ap- 
proximation apparently seems better. However, it may 
not be the case|23; higher-order terms with the ver- 
tex corrections, which are usually discarded in the self- 
consistent approximation, cancel with each other, which 
means that the lowest-order calculation such as ours is 
more reliable than the self-consistent one [2^. 

As a future problem, one should consider how to ob- 
serve the pseudogap in heavy ion collisions or proto- 
neutron stars. It is also interesting to see what would 
happen if the phase transition of the CS is strong first 
order. 
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